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The random Schrodinger equation: slowly decorrelating 
time-dependent potentials 


Yu Gu* Lenya Ryzhik* 


Abstract 

We analyze the weak-coupling limit of the random Schrodinger equation with low frequency 
initial data and a slowly decorrelating random potential. For the probing signal with a suf¬ 
ficiently long wavelength, we prove a homogenization result, that is, the properly compen¬ 
sated wave field admits a deterministic limit in the “very low” frequency regime. The limit is 
“anomalous” in the sense that the solution behaves as exp(— Dt s ) with s > 1 rather than the 
“usual” exp(— Dt) homogenized behavior when the random potential is rapidly decorrelating. 
Unlike in rapidly decorrelating potentials, as we decrease the wavelength of the probing signal, 
stochasticity appears in the asymptotic limit - there exists a critical scale depending on the 
random potential which separates the deterministic and stochastic regimes. 


1 Introduction 

We consider the weakly random Schrodinger equation 

idtfi+^Acj)-sV(t,x)cj) = 0 (1.1) 

with a low frequency initial condition </>(0, x ) = cf>o(£x). Here, e <C 1 is a small parameter measuring 
the strength of the random potential V(t, x), and i <C 1 is the ratio of the typical scale of variations 
of the potential to that of the initial data. We are interested in the long time behavior of the 
solution, on the scales such that the effect of the weak random potential is visible. When the 
temporal and spatial correlations of the random potential are decaying rapidly, this problem was 
addressed in [2] when 1 = 1 (so that the initial condition is not slowly varying), and in [T, 7] 
for t <C 1. We should also mention the papers mm where the kinetic limit for the case t = 1 is 
obtained in the much harder case of time-independent random potentials. In the aforementioned 
papers, with the rapid decay of the correlations of the random potential, the solution of (HU) is 
affected by the random potential in a non-trivial way on a universal time scale t ~ e -2 which 
depends neither on the typical scale of variations of the initial data nor on the covariance structure 
of the random potential. Moreover, the limit is deterministic for all t <C 1 - all solutions with 
slowly varying initial data homogenize at times t ~ e -2 . Here, we are interested in what happens 
when the correlations of the random potential decay slowly. As we will see, then the time scale to 
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observe “non-trivial” behaviors indeed depends on the correlations of the random potential and is 
not universal. Moreover, the observed behavior on this time scale varies dramatically depending 
on the scale of variations of the initial condition even for slowly varying initial data. 

Let us now be more specific about our assumptions on the random potential V(t,x). It is a 
stationary mean-zero Gaussian random field, with the spectral representation 



ip-x V (P, dp) 

( 2 vr) d 


and the covariance function 

R(t, x) = E {V(t + s,x + y)V(s, y)} = [ R^e^e^ 

J R d (2tt)“ 


( 1 . 2 ) 


(1.3) 


Here, 

^(P) = \ p \2j+d-2 » 0(P) = »\P \( L4 ) 

are, respectively, the spatial power spectrum of the potential and its spectral gap. The cut-off 
function a(p) > 0 is bounded and compactly supported, and /r > 0 is a constant. The role of the 
parameters 7 > 0 and (3 > 0 can be seen by setting t = 0 and x = 0, respectively, in (11.31) : 


and 


R( 0, x)= e 

J R 


ip-x 


a(p ) dp 1 

|p|27+d—2 (2vr) rf |x|2- 2 7 ’ 


|s| > 1, 


n)= f c ~^\v\ 2p t Q(P) d P 
’ ] L, \ n \2y+d-2 (O-rV 


One of the main results of (7j is that in the case 7 + /3 < 1 (rapidly decorrelating potentials), on 
the time scale t ~ e~ 2 , and after a proper phase compensation, the wave field is homogenized - 
it is nearly deterministic, as long as t <S 1. In this paper, we consider the opposite regime: the 
parameters 7 and /3 satisfy 

0 < 7 ,/ 3 <l, 7 + /? > 1, (1.5) 


so that the correlation function decays slowly. Our goal is to explore the asymptotic behavior of 
the wave field, and its dependence on the the initial condition. 

We mention that the case i = 1, that is, when the the probing signal varies on a scale comparable 
to that of the random potential, was investigated in [2]. It was shown that, after the phase 
compensation, the wave field behaves on an anomalous time scale like a lognormal distribution. 


The possible asymptotic limits 

Before presenting our main results, let us give a heuristic explanation of what one may expect. 
From now on we assume that l = e“, and investigate different values of a > 0. Consider a time 
scale e~ k with some k > 0 to be determined. The rescaled wave function 

4>S,x) = </>(A iL), 
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satisfies the rescaled Schrodinger equation 


1 _2 a -« ,• - l „,t x 


idt<j> e (t,x) + -e 2a K Acj) E — 7—j-G(—, —)cj> e (t,x) = 0, </> e (0, x) = M x )- 


The modified potential 


V £ (t,x) = - tV( —, —) 


( 1 . 6 ) 


(1.7) 


has the covariance function 

a {p) p —ii\p\ 2 ^t/e K r ip-x/e a dp _ <l(jp £'/ 1 ) r —ii\p\ 2 ^t r i(p-x)e^ a ^P 

k) ( 

( 1 . 8 ) 


R e (t,x) = 


— f 

Ik -2 / 
•7R C 


.p f-ii-i -;- = / 

e 2ft-2 y Kd |p|27+d-2 (2^)^ 7 Rd |p| 2 7+<i—2 (2ll) d ' 


which is of order 0(1) provided that the time scale exponent is chosen as 

2/3 


K = 


2/3 + 7 - 1' 


(1.9) 


It follows from m that 1 < k < 2, hence the time scale e K is much shorter than the classical 
central limit theorem time scale e -2 . We also see that for a = a c , with 


k 1 

2/3 “ 2/3 + 7 - r 


( 1 . 10 ) 


the rescaled potential V £ (t,x) has the same distribution as V(t,x) modulo the effects from the cut¬ 
off a(p). The self-similar structure of the power spectrum plays an important role in our analysis. 

Let us ignore, for the moment, the effects of the term £ 2q_k A in (11.61) . so that this equation 
formally reduces to an ODE. If a = a c , the rescaled covariance function has a limit, 


R e (t,x ) 



_fl(0)_ | p |2/3 1 ip . x _dp_ 

\ p \2y+d-2 e (2vr ) d 


( 1 . 11 ) 


as e —» 0, which can be identified as the covariance function of a generalized Gaussian random 
field W(t,x). That is, there exists a generalized fractional Gaussian random field W(t, x) so that 


E{W(t + s, X + y)W(s,y)} = Dp ^ 

Thus, when a = a c , for each x € we formally have 

d t (j) e (t,x) « —iW (f, x)4> e (t, x) 
when £ is small, so the solution is approximately 


4> e (t, x) « </>o(x) exp | — i j IT(s,a:)c/s|. 

Furthermore, when a € (0, a c ), the limit of the covariance function is 

R (t x) _- [ e -/dp| 2pt d,v 

’’ (2vr ) d J Rd |p|27+rf-2 e ap > 


( 1 . 12 ) 


(1.13) 


(1.14) 


(1.15) 
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so the spatial variable is “frozen” and there is only temporal mixing. Then, (11.61) becomes, formally: 


d t (j) e (t,x ) 

with the approximate solution 

~ -iW(t,0)c/) £ {t,x) 

(1.16) 

4> e (t,x) « 4> 0 (x) 

exp | — i 1 W(s,0)ds|. 

(1.17) 

Its Fourier transform 



M£) exp 

{ -i J W(s,0)ds} 



coincides with the limit in the case a = 0 [2] Theorem 1.2]. 


On the other hand, when a > a c , the covariance m is highly oscillatory, and so is the 
corresponding noise W(t,x/£ a ~ ac ), and it is not clear on this formal level what the limit should 
be. It turns out that the oscillations will lead to a deterministic limit. 


The main result 


We analyze the problem in the Fourier domain, where (USD can be written as 

idtA ~ ^£ 2 “ -k |£| 2 4 - V e (j) £ = 0. 

The initial profile <f>o is assumed to be of the Schwartz class: <f>o € 5(M d ). As in mm , in or¬ 
der to eliminate the large phase coming simply from the deterministic evolution, we consider the 
compensated wave function 


ipe(t,0 


4>e{t,^)e 


i\e a £\ 2 t 

2e K 


„ad 


<X- 




i\e a t;rt 


(1.18) 


Let us emphasize that £ ~ 0(1) in the argument of the function ip £ (t,^) corresponds to £ ~ 0(e“) 
in the argument of the function In the following, if we refer to the order of frequencies, it 

is with respect to the argument of 4>. 

In order to formulate the main result, we define the following constants 


roo 

K\ = Q d e~^ $ 

Jo 


dp 


where Q d is the surface area of the unit sphere in R rf , and 





| w |27+£i—2 


dw 

( 2ir) d ’ 


_ a(0)KiK 2 
~ {2n) d (2 - At)’ 

= / l s-u l ^K 2 (\s - u\t,^)dsdu. 

(2vr) .7 [o,i] 2 


We denote by N (0, a 2 ) a random variable with the normal distribution of mean zero and variance a 2 , 
and by W(t,x ) a generalized fractional Gaussian mean-zero random field with covariance () 1.12 f) . 
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Theorem 1.1. For each fixed t > 0,£ € M. d , the compensated wave function if £ (t,tf) converges 
in distribution, as e ^ 0, to a limit defined as follows: (i) if a > a c , then the limit is 

deterministic: 

fi(t,0 = 4>o{f) exp | - ^Dt 2/K Y 
(ii) for a = a c the limit is random: 

$(i, £) = [ (fo{x)e~ l ^ x exp { — i f W(s,x)ds\dx , 

J Rd 1 Jo J 

(Hi) if 0 < a < a c and (3 € (0, 1/2], or n — a c < a < a c and (3 € (1/2,1), then the limit is random: 

= <^ o (Oexp{ilV(0,T>f 2/K )}, 

(iv) if a = k — a c and (3 € (1/2,1), then 

= Mf,)^p{iN(0,D(t,f,)t 2/K )}. 


We point out that when a € (a c ,oo), the limit in Theorem II.II is deterministic, thus we have a 
convergence in probability. This means that solutions with “sufficiently low frequency” initial data 
homogenize also in the random potentials with slowly decaying correlations. The (non-standard) 
time factor f 2 / K in the exponent reflects the slowly decorrelating nature of the random media, and 
should be contrasted with the exponential decay in time exp (—Dt) of the homogenized limit in 
rapidly decorrelating potentials [7J. It is also worth noting that the wave field is attenuated in the 
limit and the exponential decay in time implies a loss of mass (the L 2 norm). The same phenomenon 
occurs in the rapidly decorrelating potentials [7], where the mass lost from the low frequencies of 
the order 0(s a ) escapes to the high frequencies of the order 0(1). 

A key feature of Theorem 11.11 is the existence of a critical wavelength scale 0(e _ “ c ) of the 
probing signal, which separates the deterministic and random regimes. This is very different from 
the Schrodinger equation with rapidly decorrelating potentials, where homogenization happens for 
all slowly varying initial data (any a > 0), as was proved in [7]. In addition, the random vari¬ 
able 1V(0, Dt 2 / K ) that arises in the limit for a < a c can be identified as the (one-point) distribution 
of a fractional Brownian motion at time t, which may be written as 

VDBi/ K (t) = [ W(s,0)ds, 

Jo 

in agreement with our informal analysis. Such fractional limits are typical for additive functionals 
of Gaussian random variables with slowly decaying correlations, but it is not a universal limit and 
comes from the specific covariance structure of the Gaussian field m- a posteriori, this limit 
justifies the basic assumption of the informal computation we have shown above: the effects of the 
Laplacian operator are suppressed via the phase compensation, the dynamics is essentially reduced 
to an ODE, and the random potential behaves as a fractional Gaussian noise in the limit. We 
should mention that the restriction a > k — a c for (3 € (1/2,1) is a limitation of the technique of 
the proof, and is the analog of the restriction f3 < 1/2 that was needed in [2j in the case a = 0 
considered there. 
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A heuristic explanation of homogenization in the “very low frequency” regime (a > a c ) may 
be as follows. Up to the time scale e~ K , the characteristic frequencies of the slowly decorrelating 
medium are of the order 0(e“ c ) - this can be seen from expression (11.81) for the rescaled covari¬ 
ance R e (t,x). For the probing signal with a bandwidth narrower than 0(e ac ), of the order 0(e a ) 
with a > a c , the interactions of the signal with the medium frequencies can only produce fre¬ 
quencies outside of the original bandwidth, so only the ballistic component survives in the limit, 
leading to the homogenization result. On the other hand, if the probing signal signal has bandwidth 
wider than 0(s ac ), that is, a < a c , then the interactions with the medium frequencies can produce 
frequencies inside the initial bandwidth, and the limit is stochastic. 

To the best of our knowledge, the first study of wave propagation in slowly decorrelating media 
was done in the one-dimensional case mm , where it was shown that a pulse going through a 
random medium with long-range correlation performs a fractional Brownian motion around its 
mean position, as opposed to the regular Brownian motion in the rapidly decorrelating case [5]. On 
the other hand, the motion of particles in such random media leading to fractional Brownian limits 
was considered in annum. We also mention the recent work of M analyzing the wave energy 
instead of phase evolution, where a time-separation is observed due to the long-range correlations 
in the random potential. 

The paper is organized as follows. In Section [2], we introduce the Duhamel expansion and 
prove some basic moment estimates. Next, we prove the homogenization result in Section [3l The 
discussion of the stochastic regimes is in Section 01 
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2 The Duhamel expansion and a moment estimate 

The basic strategy in the passage to the limit e — >• 0 is the same as for rapidly decorrelating 
potentials in [7]. We consider the Duhamel expansion for the function r ip e and study its moments 
using the expansion. In this section, we introduce the series and establish a uniform (in e € (0,1)) 
bound on the moments of the terms in the expansion. This will allow us to pass to the limit 
term-wise in the series for any moment of the compensated wave function. 


The Duhamel expansion 


A straightforward calculation shows that the compensated wave function defined in (|1.18p sat¬ 
isfies the following integral equation: 


- )ds. 

e aJ 


Mt,0 = MO + ^ f / V{ £ } g _ JL 

lE Jo J«. d V Z7r ) 

The solution of (12.11) can be written as an iterated series: 

OO 

MtiO = ^fnM’Oi 

n =0 


( 2 . 1 ) 


( 2 . 2 ) 


6 





with the individual terms 



The phase factor in (12.31) is 

n 

G n (£,s {n) ,p {n) ) = - pi - ... -Pfc-i| 2 - |£- Pi ~ ••• -Pfc| 2 )y- (2-4) 

k =1 

We used here the convention fo,e{t,Q = 4> o(£)> and have set po = 0, p^ n> = (pi,... ,p n ), as well 
as s'") = (si,..., s n ). We have also defined the time simplex 

A n (i) = {0 < s n < ... < si < t}. 


As has been shown in, for instance, mm, the series ( 12 . 21 ) converges, and one can take the 
expectation of ijj £ and its moments term-wise, as long as e > 0 is fixed. Thus, the proof of 
Theorem 11.11 boils down to the asymptotic analysis of the moments of the form 


E{/mi,e • • • fm M ,efni,e • • • fn N ,e}- 


(2.5) 


As V is a mean-zero Gaussian field, a non-zero contribution comes only from the terms with 


M N 

rij = 2k 

*=i j =i 


( 2 . 6 ) 


for some k € N. 

We will denote “the random part” in (12.51) as 


Im,N =(27 T )- 2kd V(^,dp 1A ) . . . V(^- ,d Pl>mi ). . . V(^,d P M, l) • • • V(^^-,dpM,m M ) 
e £ £ £ 


,Ul 


,UN, 1 


X dq 1,1) ... U* (^£11, dg lini )... V*(^,dq N}1 ) 


UN. 


n N 


, dq N) 


n N / 


Here, Sj and pj variables come from the terms f mk ,e, ~ each of them involves m*, temporal vari¬ 
ables S] i> i, ..., Sk,m k and momentum variables Pk,i, ■ ■ ■ ,Pk,m k , while the variables Uk and q k come 
from the terms f* k £ in (j 2.5 j) . Using the rules of computing the 2k —th joint moment of mean zero 
Gaussian random variables, we write 


E {Im,n} = J2 n e-^ w ^-^ K S( Wl + Wr )R{ Wl )^^L. 


(2.7) 


The summation extends over all pairings J- formed over the vertices 


In (12.71) . vi,v r are the two vertices of a given pair, and wi,w r are the corresponding p,q variables, 
that is, wi = Pij if vi = Sij and wi = —Qi.j if vi = Uij. The same holds for w r . We will also write a 
pair as an edge e = ( vi,v r ). Note that the order of vi,v r does not matter here since both g and R 
are even. 
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A uniform bound on the individual terms 


As we have mentioned, in order to be able to pass to the limit term-wise in the series for the 
moments of ip e we will need the following uniform bound. 


Lemma 2.1. For all e € (0,1], we have 

|lE{/ mij£ . . . fm M ,ef ni ,e • • • fn N ,e}\ ^ ^ ,, n Ar' ( 77^ 

I li=lV m V ! I 1 y = 1 \ n ]) ' 

with some constant C depending on t,(f>o,R,g, and k as in 12.6]) . 

Proof. Since cf> o is bounded, we have 

IE{/ mi , e • • • fm M ,ef* u£ • • • fn N ,e}\ < C k £^ 2k [ dsdu [ \E{I M , N }\. 

JK 2kd 

where 


( 2 . 8 ) 


W X ... X /\ nN (t ). 

By symmetry, the r.h.s. of the above expression is bounded by 

£.2(1—K,)k r r 

—jj -— „ / dsdu / |E{Im,tv}| 

n i= iK)!n ? =iK)! J[0A 2k J9. 2kd 




c k ( 2 k-iy.\ ( 2(1 _ k) 


N 


/ / 

./fo.tl 2 JS 


=-0O)b-“IA K 


[o,t ] 2 ( 27r ) c 


-R(w)dwdsdu 


V. 




with the factor {2k — 1)H coming from the total number of pairings. We recall that g(p) = /i\p\ 213 
and 

a(p ) 


and change the variable 


A(P) = 


u; i->- 


| p |2 7 +d-2 ’ 


-k/2/3 


|s — u|V 2 ^ 




to obtain 

2(1—ft) 


/ / 

Jro.il 2 ■/«' 


e - 0 (w)|s-u|/e'' 


■R(w)dwdsdu = 


e i‘l' u 'l 2/3 e K / 2 Pw dwdsdu 

—TWTm) 


l[o,t] 2 Jm d (2vr) d -v-/— 7 [o,t] 2 ,/i7 ( 2 n) d \w\ 2 ^+ d 2 “ v |s - u\ x l 2 h ’ | s _ u ^- ’ 

which is bounded since a is bounded, 7 < 1 and 7 + /? > 1, and the conclusion of Lemma 12.11 
follows. □ 

The bound in Lemma 12.11 is useful, since for fixed M,N , we have - recall, once again, that k is 
related to m* and n* via (12.61) : 


E E 


(2k — 1)!! 




mi=0 njv=0 


7‘ = E 


k =0 


(2k — l)\\C k (M + N) 2k 

J2k)\ 


< 00 


(2.9) 














as can be seen by the binomial expansion. Thus, Lemma 12.11 implies that the individual terms in 
the series for are uniformly bounded in e € (0,1] by a summable series, and 

thus not only we may write 

oo oo 

E {Mt, 0 M r e (t,Z) N }= E ••• E (2.10) 

7721=0 72jV=0 

but we may also pass to the limit individually in each summand. 


Re-writing the individual terms 

We will now look at the individual term 

E{/mi,e • • • fm M ,efni,e • ■ ■ fn N ,e\i 
and re-write in a form convenient for the analysis. We can write 

x> N 




ie* 


—IE* 


X 


r r M N (2-H) 

/ dsdu / E{I M ,N}e iQM e- iGN Y\hM,iT\ h *Nr 
JAmAt) ^ 2kd t= a j=i ’ 


where 


M N 

Qm = Y J G rnM a *i,s {mi \p (mi) )/£ K , <3 n = E G nj {e a zM ni) ,q {ni) )/e 

3 = 1 


i =1 


u Iff Pi,l + ■ ■ ■ + Pi,mi ^ u * 2*ff 9j,l + ■ ■ ■ + Qj,nj ^ 
h M ,i = -- ),n NJ = --)» 


and we recall that 


By (12.71) . we have 


A m>n (t) = A mi (t) x ... X A n N (t). 


E {Im,n} = E II e- B{wi)lvi ~ Vrl /£K 6(wi + w r )R(w t ) 

T (vi^eJ 7 


dwidwr 
(2ir) d ’ 


( 2 . 12 ) 


where T are the pairings arising in computing the joint moments of the Gaussians V. Thus, we 
can write 

E {f mi ,e ■ ■ ■ fm M ,ef ni ,e • ■ • fn N ,e} = E J m 1 ,...,nj v (^), (2-13) 

T 

with the individual terms 


T £ 

u m 1 ,...,n* N 


(D= 


e \E“,m f £ 


is* 


x / dsdu / 

JA m , n (t) J^ 2kd /- v 


M N 


o-B(wi)\vi-v r \/e r 


5(wi + w r )R(wi)e lGM e iGn h M ,i h 


dwidw r 


N,j- 


(vi,Vr)&F 


i= 1 j =1 


(2vr) c 


(2.14) 
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Our goal is, therefore to identify 


limJl n * (J 7 ), 

e ^0 lr "’ JV 


for a fixed pairing T. We integrate out the w r variables in (I2.14p and obtain 

x—^ N 


= 7 


IE n 


— l£ n 


M N 


x[ dsdu [ TT e- B ^\ vl - Vr \/ £K R(wi) e iGM e - i6N f\h M ,if\h^ 

JAm,n(t) JR kd ( VuVr ) e jr i=l j=l 


dwi 

J (2ir) d ’ 


with Qm-, Gn, h,M,i, hwj subjected to the constraint 

wi + w r = 0 


(2.15) 


(2.16) 


for all pairs. Once again, we recall that 

0 (p) = mM 2/3 , R(p) = j p | 2 7 +L 2 - 

Using these expressions, for every wi, we apply a key rescaling to extract the characteristic frequen¬ 
cies from the underlying media up to a time scale e~ K \ 

Wl !->• £ ac Wl. 

Using the fact that a c = k/2/3 and k = 2(3/(213 + 7 — 1), we obtain 


T £ 

u m 


(, J-) = - - - - - 

'f l Z-^i=l rn i ( _ 7=1 


3 = 1 ,L 3 


x / dsdu / TT e 

JjR kd ( VltVr ) e j: 


( a \ M JV , ( 2 -17) 

^ Ell p . p 


| Wz |2 7 +d-2 


*=1 J=1 


where 


M JV 

<5jvf = ^2 G m i ^ a ts (mi \e a ^ rn G)/ £ ^g N = J2G nj (£ a Z,u<- ni) ,£ aB q in * ) )/e K , 

l=i 


(2.18) 


Z=1 


and 


hM,i = MC ~ £ ac a (Pi, 1 + ■ ■ ■ +Pi,m;)), h* N j = 4> o(£ - e“ c “(<?!,1 + • ■ ■ + Qj , nj )), (2.19) 

which are still subjected to the constraints (12.161) for all edges. This expression will be the starting 
point for our analysis for all frequencies. 
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3 Homogenization of the very low frequencies 

We begin with the very low frequency regime (ct > a c ), where the limit does not depend on 
whether j3 < 1/2 or /3 > 1/2, and prove the homogenization result - the limit of the compensated 
wave function is deterministic. 

Proposition 3.1. If a > a c , then, for a fixed t > 0 and £ G M. d , we have 

-> M0e~* Dt2/K 

in probability as e —>■ 0. 

Proof. We will first compute the limit of E{?/> e (t, £)}, and then show that the second absolute 
moment converges to the square of the limit of the first moment, which will establish the convergence 
in probability. 


The limit of the first moment 


Let us start with E{^ e (t, £)}. In that case, there is only one simplex of time variables, and (12.171) 
simplifies to 


J £ 

J n 


(F) = 


4 MO [ ds f JJ 
1 Ja 2k (t) JR kd /T. 


a(£ ac Wl) iGn(e a i,s^ n \e ac v^ n ' > )/e K dvjj 


(vi,Vr)&F 


\wi\ 2 ^ +d 


-2 


(2vr) 


with n = 2k for some k G N. By the definition (12.41) of G n , and since a > a c = k/ 2/3, we have a 
bound for the phase factor 

|G n (e"e,s (n) ,e“ c P (Tl) )/e K | < Ce%~ K , 

with a constant C depending on £, s( n ),p( n ). Since f3 < 1, the function a(p) is uniformly bounded 
and the function 

e -a\v\ 2 p \s-u\ /|p| 2 'v + rf — 2 5 

is in L 1 ([0,t] 2 x M rf ), we may apply the Lebesgue dominated convergence theorem. This gives 


JkiJ 7 ) 


a 2 k 


MO ! ds [ TT °(°) 


Therefore, we have 


lirn E{ip £ (t, 01 = 5^ !™n ^fcG 77 ) 


/c=0 J 7 
oo 


£—^0 


/» 

££(-d‘a>(o 

i_ r\ T~ J l— 


k =0 J~ 


ds JJ 

a 2 ,W ^ M (w , )tv)GF 


3 —tv | ^(0) _ dwi 

\ Wl \2n,+dr-2 ( 27r )d- 


Integrating out the wy variables, and recalling the definition of 


K x 



dp 


ll 









we obtain 


liniE {^ e (t,0} = 

k =0 


a(0)Ki 

(2vr) c 


i fc 


E 



— f r 


i—y 

P . 


Note that 

X] II I vi-v r \~~^ 

F {vi,v r )£T 

is symmetric in (si,..., S 2 k), so 


x / ds n \ vi ~ Vr \ pl 

T ^2k(t) 


(2k-l)\\ 



i —i 

s— u\ P dsdu 


and thus 

limE{V> e (t,0} = ^o(O e_ ^ Di2/ "- 


The limit of the second absolute moment 

Next, we consider E{|'0 e (i,£)| 2 }, then ()2.17j) becomes 

dsdu / 


r m =J_X_ f 

) ■rr > / ■\ ri I 

i m {-i) n j A 


A m,n (^) 


n « 

(vi,V r )&T 


-H,\ m \>P\v r -v r \ a(£ ac Wl) dwi 
\wi\ 2 ^ +d ~ 2 (27r) rf 


xe iGrn{e a Z,s( m \£ a c p {™))/ £ K e -iG n (£ a Z,u( n ),e a c q ^))/ £ K _ £ a c -a ^ p .)<£*(£ 


2—1 


- £ “'"“E«j)- 

i=i 


with m + n = 2k for some k € N. Here, the p, g variables satisfy the constraint 


u>z + tty = 0 


for all edges. 

The “crossing” pairings that have at least one edge e = (vi , v r ) such that vi is an s- variable, 
and v r is a u- variable satisfy 

m n 

X* = -X^°- 

*=i i=i 

Thus, such pairings give the following factor coming from the initial condition 

m n 

0o(C-e ac_ “X K ^°^ -e “ C_a X] <7 ^ 0 

i =1 J=1 

as e —)• 0 , as a. > a c . If we denote the set of those “crossing” pairings by P Cl then it is clear that 


lim 

£—^0 


X <4,n*cn = o. 


TdPc 
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For T £ P Cl all edges e = (■ vi,v r ) connect the same type of variables, thus both m and n are 
even, and 

m n 

Y P i = Yv = 0 ’ 

i= 1 j =1 

under the constraint wi + w r = 0. Then, by the same argument as for E{^ e (t,£)}, we have 




1 1 


i m (—i) r 


(0 [ dsdu [ 

J dRfcd , 7 


1 1 

"i™ iPi) 


777-|0o| 2 (O 


A m,n (t) 

a{0)Ki ik 
. (2ir) d 


(vi,V r)€-F 


= —u r | a (Q) _ dwi 

\wi\^+ d ~ 2 (2vr) d 


/ dsdu n h 


1~7 


— v r 


>Am,n(t) {vi,V r )&T 

Since m, n are both even, we write m = 2ki,n = 2/^ for fci, &2 £ N, and obtain 


limE{|V> e (t,£)| 2 } = ^fci^O 77 ) 

^1 5 ^ 2—0 TtflPc 


■ E (-i)‘ ,+ ‘ 2 i*)I 2 (q 

ki,k2=0 


a( 0 )Ki i ^1+^2 


(27) 


s n h - u r | 

T$pJ Am - r *( t ) 


1~7 


As P P c , we have 


1-0 


/ dsdu JJ |u; — u r | 0 

jr^pJ A mt n{t) {v U Vr)&T 


y ds n ^ 

■A (v u v T )£F a 


, 1~7 

— V J 0 lx 


£/ a n i”‘ 

■Ax 2fe 2 ( ) {yi,V r )&Tu 


, 1-K 

— uJ 0 


where P S ,P U denote the pairings formed by si,..., S 2 fci and ui,..., U2fc 2 respectively. This implies 

hmE{|'0 e (t,£)| 2 } = I limE{^> e (t,£)}| 2 , 

e-fO £->-0 


which completes the proof. □ 


4 The stochastic limits 


In this section, we assume a € (0, a c ] and prove a convergence in law of the compensated wave 
function ijj e (t,£). To identify the limiting distribution, we compute the limiting moments. Recall 
that by Lemma 12.II we only need to pass to the limit in (12.171) : 




(P) = 


E m _ 

i=i m i ^_j j yZj=i n j 


/ dsdu / 

JR kd , 


C Vi,Vr)eT 


,-^\ Vl -v r \ <^l) dwt p iQM e -iQ N jy hM . -Q ~ h * N j ' (4 ‘ 1} 

i= i i=i 


|^|2 7 +d- 2 ( 2 7r) d 
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The phase factors are 


M 


N 


Qm = Y J G mM a ^s^\e a ^)/e\ g N = Y J G ni {e a ^u^\e^q^)/e\ 


i =1 


thus for a < a c , we have 


3 = 1 

\Qm\ + \Qn\ ~ e a+CXc ~ K . 


We will need to consider two cases. 

Case 1: /3 € (0,1/2]. Then a c > k, thus we have 

a + a c — k > 0 , 

which implies Qm- Qn —>• 0 as e —>• 0, so by m 

i i 


r 

U ryy 




E m 

i=l m i = 


(-0- 


-j=l “J 


[ dsdu [ n e-Mkzl^hz-.d °(°) * wi H a &p,q), 


(4.2) 


with 


M N 

H a (£,P, q ) = lim h M ,i 


(4.3) 


»=i i=i 


J (<A)(OrWo(0)" ^ a € (0, a c ), 

\ n"i 0o(£ - Pi,i - ... - Pi, mi ) nf=i kit - Qj, 1 - • • • - Qj,nj) a = a c . 

We recall that all p. q —variables satisfy the constraints wi + w r = 0 for all edges. 

Case 2: /? G (1/2,1). Here, we will consider the range 


k — a r < a < a r 


We have 


or 


a + a c — k > 0 for a € (k — a c , a c \, 


a + a c — K = 0iora = K — a c . 
When a € (k — a c , a c ], we still have 

Gm iQn —^ 0 as £ —> 0, 

so we have the same limit as (14.21) : 

1 1 


T £ 


(H) -i 


E m _ ^ 

i=l m i ^_j j ylj=l n 3 


[ dsdu [ TT e -Awi\ 2 f> \vi-v r \ ^(0) _ dw I . , 

L mAt) U Jl/ N 2 ^- 2 (2- y Ha{ ^ q) - 


(4.4) 
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On the other hand, when a = k — a c , we have 


M mi N Tij 

Qm ->■ ^ ^{Pi,k • 0 s i,fc and ->• • O u j,k, 

i =1 fc= 1 j=l fc=l 


hence 


T e 




E m ^ \-~^N 

i=i rn i ^i 


- / (Is (I'll / 

! n J JA m , n (t) JwL kd "7. 


(vi,v r )eT 


^ u \ wl \ 2 P\ vl _v r \ g(0) dwi 

\ W l\2'y+d-2 (27T ) d 


To summarize, depending on the parameters a, /?, 7 (/c depends on ,5,7), we can write 


(4.5) 


hmJ^.^GF) — 

with </ mi ,.... n * v (F) given by (14.21) . (|4.4|) and (14.51) in the three respective cases: (i) 0 < (3 < 1/2 
and 0 < a < a c , (ii) 1/2 < /3 < 1 and k — a c < a < a c ; and (iii) 1/2 < /3 < 1 and a = k — a c . 


The case a < a c 


Theorem o treats two regimes in the case 0 < a < a c : (i) f3 € (0,1/2] and a € (0, a c ), and 
(ii) /3 € (1/2,1), and a € [re — a c ,ot c )- The key feature in this case is that the factor (14.31) coming 
from the initial conditions in (14.ip has the limit 

M N 

*=1 3 = 1 


as e —>• 0. If, in addition, we assume that either /3 € (0,1/2] (case (i) above), or f3 € (1/2,1) 
but a 7^ k — a c in case (ii), then the term J also simplifies: the phase information vanishes in the 
limit, and 


J, 


mi 


A?) = 


1 


E m x ~^n 

i=l m i = l n j 


/ dsdu / 

J A m ,n(£) J K 


(4.6) 


, n _ 2 77 (Monmr 

{vi,V r )&T ' 11 K ’ 


Making a change of variable wi 1—?> wi\vi — v r \ 1 ^ 2/3 and integrating out wi , we obtain 


Jm\,...,n* N (-E) — 


jEi=i m i (_^E/=i%- 

a(0)ivi i k 


(2tt) c 


(Monmr [ dsdu n 

(vi,v t )&T 


— Vr 


1-7 

P . 
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As in the proof of Proposition 13.11 we use symmetry after summing over all pairings J- to get 


y ] Jmi,...,n* N (J~) 


1 


1 


T 


which implies 


a(0)Ai 

*E£ H) Ef = i^ L (27T Y 

x (r 

nZimmf =1 nj\ vw 


oo oo 


l k 




. _ 1~7 

s — u\ P dsdu 


lhnE {Mt,0 M Tp*At,Z) N } = - 

^ n iEi=i"u f_ ? -'|Ej=i%- L (2vrJ a 
m i=0 n N =0 l \ l J 


(M( » M (m)) K 


X 


(2k - 1)!! 




1-7 \ 

|s — u\ P dsdu , 


M 


N 


V 


with 

2A: = ^2 mi + n, 
i =i l=i 

in the summand. The binomial expansion tells us that 

^ 1 1 (2fc)! ,M N, 

,-Eili "U ( ^Eli n, n M ., 

mi+...+njv= 


=2fc *Ei= imi (-i)EiLi «* Hill mi! ni=i 


= (— - — ) 2k = (M - NY k (-l) k , 


thus we have 


limE {Mt,S) M r e (t,Z) N } = J2(M - N) 2k (-lY 


o( 0 )/m 


i fc 


k=0 


L (2tt) c 


1 


(*>(0) M (to* 

i 


2 k k\ 
= E 


'[o,f 

.(Oe 


| s _ nl-Vdsdu = (^o(0) M (^(0) iV exp { - -(M - iV) 2 Ht 2 / K } 


iN(0,Dt 2 / K ) 


M 


~iN(0,Dt 2 / 


■>n 


which completes the proof in the cases 0 £ (0,1/2], a. £ (0, a c ) and f3 £ (1/2,1), a £ (k — a c , a c ). 
In the case 0 £ (1/2,1 ),a = k — a c , the phase does not disappear in the limit, and we have 

1 1 


Jm\,...,n* N (d~) — 


E m r _ v "'n 

i =1 m i ^_^E7 


- / dsdu / 

i"i iA m ,„(t) J'R kd t_. "7 \. 


( Vl,V r )£J r 


„-u\wi.\ 2 P\V1-Vr\ Q (°) dW l 

|^|2 7 +d-2 ( 27r )d 


x e * E"i ESi(w.fc-0*i,fc e -< Ef=i Efcii(©,*-f)«,■.* (0 O (<£)) M (0o (£)) ^ 

Compared to (|4.6I) . the only difference is the extra phase factor 

M mi N n j 

j=l fc=l 1=1 fc=l 


(4.7) 
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Since the p, q —variables in the above expression satisfy the constraint 


wi + w r = 0, 


we can write (El) as 
1 


J, 


m\ 


XJ 7 ) 


iy [ dsdu [ n e 

1 d A ^ Rkd (vi,V r )eT 


dvn 

\ W l\2 1+ d-2 (2vr) d 


(-*) 

x (MC)) M (m)) N - 

Once again, we may change the variable wi i->- wi\vi — v r \~ 1 / 2 ' 3 and integrate out wi to obtain 

1 1 


J, 




AD = 


Kor(0o«)) M («{))' v 


E m _ v ™\N 

i=l m i 

/ dsdu JJ |u/-u r | ~AT 2 (h -u r |,£)- 

J \ s- T~ 


We recall that 


#2(A,0= [ e 
J R 


^ (vi,v r )eT 

..( io - OA 1 -^) -1 dw 


—li\w\ 2 P 


| w | 27 +d 2 (27r) rf 

By symmetry we have, after summing over all pairings 

1 1 


(4.8) 




[a(0)] fe (<^o(0) M (^(0) JV 


J 7 


E m v—viv ^ 

i=i m * ( — 7*^J=l n d 


x 


(2fc- 1)!! / r 

nf =1 ^!nL^! [MV 


1~7 


s — u| 0 dC 2 (|s — u|, £)dsdu , 


\7V 


and as before, this implies 

OO 

lim E {Mt,() M rAt,() N } = E< M - A r ) 2 ‘(-l) l [a(0)]‘(*K))"(*(C)' 

£—>-0 L ' 

k =0 

x ¥v. (/ otp i»-= (uo) M (m)) N ^ M - N?DM,v ' 

The proof in the case (5 € (1/2,1), a = k — a c is now complete. 

The case a = a c 

When a = ct c , the argument in the initial condition in the expression for J is different: 




1 


1 


E m 

i=l m i (_ 

M 

-[M€ -Pi,i 


N - / dsdu / 

pi"! JA m ,„(f) 


(yi,v r )&F 


^-p.\wi\ 2 P\vi-y r \ a (°) dw; 

|^|2 7 +d-2 ( 2 7T) d 


AT 


X 


i=l 


Pi,mX tJ '^o(‘? 
i=1 




(4.9) 
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The goal is to show that 


oo oo 


where 


limE mt,o M rAt,o N }= E ••• E 

mi=0 njv=0 T 


il>(t,£)= [ (po(x)e l ^' x exp { — i f W(s,x)ds\dx, 

jRd 1 Jo J 


and W(s,x) is a generalized Gaussian random field with the covariance 


E {W ( s , x)W ( u , y )} = lZ{s — u,x — y) = e 

J ffi' 

Let us define 


—/i|io| 2 ^|s— u\ iw-(x-y) 

| yj1 2 7+d—2 (2n) d ' 


fn(t,0 = - 7 / 00 {x)(—iVt{x)) n e ^' X dx 
n\ Jud 


with 


Vt(x) = f W(s,x)ds. 

Jo 

The following lemma will help us express the moments of ij: as a series expansion. 

Lemma 4.1. We have 

OO 

and 


(4.10) 

(4.11) 


n=0 


oo oo 


e mo M nt,o N }= E ••• E 

mi=0 npf=0 


(4.12) 


Proof. For every t > 0 fixed, Vt(x) is a mean-zero stationary spatial Gaussian random field, so 
for fixed t > 0 and x € we have 


E -(-iV t (x)) n ^ e~ iV ^ 

,t^o re! 

in L P (Q) as N —>• oo with the L p error independent of x for any p > 1. This implies 

N 

^ ~2fn{t ,£) -t 

n=0 

in L p ( fl) as N —> oo, which further leads to the moment representation (14.121) . □ 
For fixed mi,..., njv G N, assuming that 

M n 

E m * + E rij = 2k, 

i =1 j=l 
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for some k € N, we write 






(4.13) 


n,! 


1 't? 


Since 


we have 


„ M N 

7 R (M+JV) d AJ. 

r. m „ m 

V t (x) m = / TT W(si,x)ds = ml / TT TU(sj,x)<is, 




E M v^lv 

i=l m i (_j)5Dj = l n j 


r M N 

/ dxdy TT 0 O (xi)e^' Xi TT 0o(y 

y M (M +J v )d ^ y 


x E 


M mi N n j 

n n w(sifa, x i)d s i tkl jj n 

Am,n(t) i=1 fcl=1 


=>*£■% 


(4.14) 


1=1 *2=1 


Using the rules of computing the 2k —th joint moment of mean-zero Gaussian random variables, we 
obtain 


M mi 


N rij 


e (q n w(siM> x i) 

i= 1fci=l 


n n ^( u i,fen%-)} 

1=1 fe 2 =i 


n n ^ ^ ^=z) i hd n 

•T 7 T 


-fi\wi\ 2p \vi-v r \ a (0) iwy(zi-Zr) dwi 

\ Wl \2 1+ d-2 c (2vr) rf ’ 


C Vl,Vr)€T 

where the summation extends over all pairings T formed over the vertices 

{^ii, . . . , Si i?7ll , • • • , SM,1, • • • > SM,mM^ 1,1 > • • • i Hl,ni j • • ■ ? HJV, 1 j ■ • ■ j }• 

Here, v/,v r are the two vertices of a given pair, and zi,z r are the corresponding x,y variables. 
Now, (14.14|) becomes 


E {fm 1 ...fm M f: i ---fn N } = 

M N 


1 


1 


E m v - ' iv „ 

i= I m * (— iyli=i n i 


y / dsdu / dxdy 

j JA m , n (t) J R (M+N+k)d 


X 


n«*i)e- i5 « n«(w)e i5 « n e -" 1 - 1 


*=1 1 = 1 (vi ,V r )£T 

Since each wi corresponds to a given pair (vi,v r ), we introduce rc r and write 


(2tt) 


under the constraint 


e iwi-(zi-z r ) _ p iwizi e iw r z r 


wi + w r = 0. 


(4.16) 
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We say that wi corresponds to vi and w r corresponds to v r . If a w variable corresponds to some Sij, 
we denote it as pij ; if it corresponds to some Uij, we denote it as —qij- Therefore, we can write 


JJ e iwiZl e iw r Zr 
{vi,v r )tT 


n 


N 


e »®*-(p»,l+—■ +Pi,m 





2—1 


3 = 1 


and an integration in x, y in (14.151) leads to 


■ ■ ■ fm M fni • • ■ fn N } 

1 1 


E m _ v—ov „ _ 

i=l m i (—^ jA m>n (t) 

M N 

x n <^o(^ - Pi,i - • • • - Pi.mj n - 9j,i 


E 


dsdu 


n < 

(vi,Vr)£F 


—[i\wi\ 2 P\vi—v r \ a (Q) _ dwi 

\ Wl \2r,+d-2 {2 -kY 


(4.17) 




i =1 j=l 

under the constraint (|4.16l) for all (vi,v r ). Comparing to (|4.9|> . we see that the proof is complete. 
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